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1. INTRODUCTION 

Bose-Einstein condensation is the macroscopic occupation of a single quantum 
state. In bulk liquid Helium, most theoretical calculations [1] agree that, in the limit 
of zero temperature, roughly 10% of the Helium atoms are condensed in the ground 
state. For recently observed atomic condensates in 87 Rb [2], 7 Li [3] and 23 Na [4], 
nearly 100% condensation is expected because of their low density. This collective 
occupation essentially magnifies ground state properties by a factor N equal to the 
number of condensed atoms. For macroscopic iV « 10 24 , this exceedingly large factor 
can be exploited to detect minute changes in the condensate ground state. Thus aside 
from its intrinsic interest, BEC is potentially a very sensitive probe. 

In this work, we will explore changes in the condensate ground state induced by 
a "sizable" impurity. By "sizable" , we mean an impurity whose size is within a few 
orders of magnitude of the trap size, and not necessarily atomic in scale. A sizable 
impurity will "drill" a hole in the condensate wave function and alter its energy. The 
question is whether this microscopic change can be detected macroscopically because 
of the Bose-Einstein condensation effect. This is described in the next section. We 
discuss the effects of interaction in Section 3 and the case of two or more impurities 
in Section 4. 



2. IMPURITY EXPULSION 



The effect we seek to describe is generic. We will therefore consider the problem 
in its simplest conceptual context. We imagine a non-interacting Bose gas of mass 
m confined to a spherical cavity (an infinite square well) of radius b. We introduce 
an impurity, which to first order approximation, can be regarded as a hard sphere 
of radius a. When the impurity is placed inside the cavity, the wave function of the 
Bose gas must vanish on the surface of both the cavity and the hard sphere. Since the 
cavity only serves to confine the Bose gas, it may or may not also trap the impurity. 




Figure 1. The Bose gas is confined in a spherical cavity of radius b and 
excluded from an off-center impurity, which is a hard sphere of radius a. 



Classically, since there is no interaction between the cavity and the hard sphere, the 
hard sphere is free to be anywhere inside the cavity. However, when the cavity is 
filled by a Bose gas, the position of the impurity is dictated by the ground state 
energy of the Bose gas. Thus there is a quantum mechanically induced interaction 
between the hard sphere and the cavity. 

The situation is as shown in Fig. 1, where the hard sphere is displaced by a 
distance d from the cavity center. The impurity is assumed to be sufficiently massive 
that the Born-Oppenheimer approximation is adequate. The effective potential for 
the impurity is then just V(d) = NEo(d), where Eo(d) is the ground state energy of 
a single particle confined in a spherical cavity with an off-center hole: 

'* -V 2 t/)(r)=-B(#(r), ip(r) = at |r| = a and |r + d|=6. (1) 



2m 

When d = 0, both the ground state wave function and the energy are known analyt- 
ically 

Mr) = - sin -rr ^ , (2) 
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(b-a) 



E °~2^Jb^W (3) 

When the impurity is off-center, spherical symmetry is broken and the problem is 
non-trivial. Instead of solving the problem exactly, we gauge this off-center effect by 
a variational calculation using the following trial function 

0(r)=(l-^)(6-|r + d|). (4) 



The first factor is the zero energy solution (Laplace's equation) satisfying the Dirichlet 
condition at Irl = a. The second factor forces the wave function to vanish at the 
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Figure 2. Left: The ground state energy of a particle in a spherical cavity of 
radius b with a displaced hard sphere of radius a. The energy as a function 
of the displacement for various impurity sizes is normalized by dividing by the 
exact on-center energy. Right: Similar calculations for the case of an impurity 
in a harmonic well. In this case b is the harmonic length and the energy is only 
normalized by dividing by the variational on-center energy. In both cases, the 
impurity is not assumed to be trapped by the cavity or the well. 



displaced cavity surface. For d = 0, the resulting variational energy simply replaces 
the factor n 2 in (3) by 10, which is an excellent approximation. For finite d, the 
resulting energies are as shown in Fig. 2. The angular integration is sufficiently 
cumbersome that we computed the variational energy by the method of Monte Carlo 
sampling. The energy is lower as the impurity moves off center. This result is easy 
to understand once one has seen it. The impurity "drills" a hole in the cavity's 
ground state wave function. It is more costly to drill a hole near the wave function's 
maximum (the center) than at its minimum (the edge) . Thus a sizable impurity will 
be expelled by a Bose condensate. 

Similar results are obtained if one replaces the spherical cavity by a spherical 
harmonic oscillator. The size paramenter b can be identified as the harmonic length 
via 

h 2 l 

hu = (5) 
m b z 

The harmonic oscillator ground state wave function is then 



^o(r)=exp[-^) 2 ]. 



(6) 



When an impurity is placed in this harmonic well, the ground state wave function is 
again altered. We now do a variational calculation with the trial function 



^(r) = (l-^)exp[-i(!l±^], (7) 

where a is a variational parameter. The angular integral can be done exactly when 
we displace the harmonic oscillator rather than the impurity. The remaining radial 
integral can then be computed by numerical quadrature. The results are also shown 
in Fig. 2. In order to compare this with the cavity case, we must take at least 26, 
and not just 6, as the equivalent cavity radius. Since the exact on-center energy in 
this case has no simple analytic form, we normalize by dividing by the variational 
on-center energy. 

To get an order-of-magnitude estimate of the expulsion force, consider the case 
of lithium atoms with ^- k 7 K A 2 and b ~ 3 x 10 4 A. The characteristic energy 
scale is 

Ti 2 1 7 

?iOJ = — -r W j-tt « 10" 8 K. 

m b 2 (3 x 10 4 ) 2 

The variation in energy is over some fraction of the trap radius. Thus the expulsion 
force due to each atom is 

„ 10" 8 K 10" 8 K10- 13 J/K 2 . T 

F « j-s- « — —o — « 10" 25 Newton. 

10 4 A 10 4 A10- I0 m/A 

For N « 10 24 Bose condensed atoms, the expulsion force IS macroscopic. Currently, 
the maximum N achieved in any of the atomic trap experiment is only about 10 6 . 
Such a small force may still be detectable if the impurity is sufficiently light. 



3. THE EFFECT OF INTERACTION 

To gauge the effect of interaction among the Bose condensed atoms, we do a 
variational Hartree calculation. For interacting Bosons confined in a harmonic well, 
the Hamiltonian is 

* = £b-v? + £>+£»(*«>■ w 

i i>j 

In the low density limit, the two-body potenial can be replaced by the scattering 
length approximation 

h 2 

virij) ^4tt— a sc 5 3 ( riJ ). (9) 
m 

For a Hartree wave function consisting of a product of normalized single particle 
states <j>(ri) 

n 

*(r 1 ,r 2 ...r n ) = J]0(r i ), (10) 
i=i 

the variational energy is given by 

W = £ [/ rf3 ^( r )(" V2 + + ( N ~ i) 4 ™- / ^ 2 (r)^ 2 (r)] . (11) 
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Figure 3. Left: The change in energy at one impurity size a = 6/5 as a func- 
tion of the interaction strength g = Na sc /b. Right: The change in energy for 
various impurity sizes at one negative interaction strength g = —0.65. Varia- 
tionally, the system would collapse at g = —0.67 without any impurity. The 
straight segments correspond to impurity locations that would destabilize the 
condensate. 



Minimizing this with respect to 0(r) yields the Gross-Pitaevskii equation. To account 
for the hard sphere impurity, one must again require 0(r) to vanish on the impurity's 
surface. Instead of solving this problem exactly, we simply take 

^ = 7^ {1 -¥TJl )eM -^ n (12) 

where Z is the normalization integral, and minimize the energy functional(ll) with 
respect to the parameter a. This is in the spirit of using Gaussian trial wave functions 
to study the Gross-Pitaevskii equation, as suggested by Baym and Pethick [5] in the 
context of BEC. 

The interaction is characterized by a strength parameter g = Na sc /b. The left 
panel of Fig. 3 shows the effect for positive scattering length. For g < 5, the results 
are only slightly higher than those in the last section. All atomic experiments are far 
below the g = 5 limit. In the extreme case of g > 10, the picture changes completely 
and the impurity is stabilized at the center. However, when the effective interaction 
is this strong, the mean-field approximation is no longer creditable and one must 
consider two-body correlations, as in the case of liquid Helium. 

The effect of negative scattering length is shown on the right panel of Fig. 3. 
Variationally, as shown by Fetter [6], the condensate by itself is unstable below 
g < —0.67 . When g is close but above this critical value, the introduction of a 



sufficiently large impurity will cause the the condensate to collapse. This is shown 
by the straight line segements on Fig. 3. For these locations of the impurity, there 
are no energy minima for the variational parameter a. When the impurity is at the 
center, it increases the single particle energy and there is no collapse. However, as 
it moves off center, the energy is lowered and the condensate collapses. This is also 
understandable from the opposite direction; when g is close to the critical value, as 
the impurity enters the condensate, it increases the condensate density without sub- 
stantially increasing the single particle energy. It pushes g over the critical value and 
the condensate collapses. 



4. TWO OR MORE IMPURITIES 

For the case of n impurities, each having a different radius ai and located at d^, 
one can simply generalize the single particle trial function (12) to 

1 n 1 

^ r ) = 7i n^w ex p[-2 ( ^ )2 ]' (i3) 

^ i=l 

with 

The resulting expression for the variational energy has a simple form 

-2 



Ey _ %' 
N m 



where 



i=l i=l 



(16) 



and the expectation value is with respect to the trial function (13). The first two 
terms in the expectation value involving may be interpreted as the kinetic energy 
of each impurity and their collective interaction with the harmonic well. With more 
than one impurity, there is no way to do the angular integration exactly. We evaluated 
the expectation value in (15) by the Monte Carlo method. 

For n = 2, the resulting energy E2 = Ey/Tiuj is shown in Fig. 4. In order to 
disentangle the two-impurity interaction energy from the effect of off-center energy 
dependence, we separate the two impurities by keeping each at equal distance from 
the trap center. The configuration used is as shown on the left panel of Fig. 4. The 
resulting energy on the right panel clearly shows that there is an effective attraction 
between the two impurities. This is induced by the Bose condensate. It is less costly 
to "drill" a slightly larger hole at one place than to "drill" two separated holes. 
One can therefore infer that multiple impurities tend to clump together and will be 
expelled from the trap center. 

The effect of interaction can again be assessed by including a Gross-Pitaevskii 
type mean-field interaction. For the present qualitative discussion, we have not both- 
ered to include this correction. 
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Figure 4. Left: Two impurities maintaining the same distance from the har- 
monic well center while rotating through an angular separation 9. Right: The 
energy for two impurities, each of size a = 0.26, as a function of angular sepa- 
ration at various distances from the well center. 



5. CONCLUSIONS 

In this work, we have considered the possible use of BEC as a sensitive probe 
for detecting microscopic changes in the condensate ground state. Our variational 
studies suggest that 

a) A hard-sphere-like impurity will be expelled from the center of a condensate. 
For a light but sizable impurity, such an expulsion maybe macroscopically 
observable. 

b) A sizable, hard-sphere-like impurity will accelerate the collapse of a conden- 
sate with negative scattering length. 

c) A Bose condensate induces an effective attraction among hard-sphere-like 
impurities. This may have interesting implications for induced dimerization 
or clusterization of weakly interacting impurities, such as 3 He. 

Work is currently in progress to seek exact solutions for problems that have only 
been solved variationally in this work. 
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